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ANALYTICAL  STUDY  OF  DIFFUSION -CONTROLLED  CHAR  OXIDATION 


AND  ITS  EFFECT  ON  STEADY-STATE  ABLATION 
OF  PLASTIC  MATERIALS 

By  Gerald  D.  Walberg 
Langley  Research  Center 

SUMMARY 

(An  analysis  of  diffusion-controlled  char  oxidation  has  been  carried  out  and  differs 
from  previous  analyses  in  one  important  respect:  previous  investigators  have  assumed 
that  only  the  carbon  present  in  the  char  would  react,  but  the  present  analysis  has  shown 
that  the  carbon  present  in  both  the  char  and  the  pyrolysis  gases  (i.e.,  all  the  carbon  orig¬ 
inally  present  in  the  virgin  plastic)  reacts  and  that  the  reaction  rate  is  determined  by 
the  rate  at  which  oxygen  can  be  supplied  by  diffusion  through  the  boundary  layer  plus 
pyrolysis  of  the  plastic.  Because  of  this  fact,  the  char  oxidation  rate,  the  steady-state 
ablation  rate,  the  boundary -layer  composition  at  the  char  surface,  and  the  heat  released 
by  chemical  reactions  were  found  to  be  dependent  on  the  composition  of  the  virgin  plastic 
and  the  pyrolysis  gases^ 

For  a  given  flow  environment,  the  char  oxidation  rate  and  the  steady-state  pyrolysis 
rate  were  found  to  increase  as  the  hydrogen,  nitrogen,  and  oxygen  contents  of  the  virgin 
plastic  increased  and  to  decrease  as  the  carbon  content  of  the  virgin  plastic  increased. 

It  was  found  that  the  heat  released  by  both  gas  phase  and  heterogeneous  reactions  can  be 
accounted  for  by  an  "effective  enthalpy  potential"  which  is  a  function  of  the  stream  oxygen 
mass  fraction  and  the  compositions  of  the  virgin  plastic  and  the  pyrolysis  gases.  Partic¬ 
ularly  high  heat  releases  are  predicted  for  polymers  which  contain  nearly  equal  numbers 
of  carbon  and  oxygen  atoms  and  produce  pyrolysis  gases  containing  large  quantities  of 
C2H2  and  C2H4. 


INTRODUCTION 

Since  the  efficiency  of  a  charring  ablator  is  determined  to  a  large  extent  by  the  char 
removal  rate,  an  understanding  of  char  removal  mechanisms  is  essential  for  the  predic¬ 
tion  of  the  performance  of  a  charring  ablator.  For  the  steady-state  ablation  case,  the 
linear  rate  of  pyrolysis  is  equal  to  the  char  removal  rate;  hence,  the  virgin  plastic 
pyrolysis  rate  is  known  if  the  char  removal  rate  can  be  determined. 


When  the  char  layer  under  consideration  possesses  adequate  structural  strength  so 
that  mechanical  failure  (removal  by  aerodynamic  shear,  spallation,  etc.)  does  not  occur, 
the  char  removal  rate  is  determined  by  the  rate  of  chemical  reaction  between  the  char, 
the  gaseous  pyrolysis  products,  and  the  chemical  species  present  in  the  boundary  layer. 
Since  an  ablator  which  loses  large  amounts  of  char  due  to  mechanical  failure  is  not  suit¬ 
able  for  use  on  reentry  heat  shields,  it  is  not  unreasonable  to  expect  that  chemical 
reactions  will  be  the  primary  char  removal  mechanism  during  most  reentries. 

The  maximum  rate  at  which  char  oxidation  can  proceed  is  ultimately  determined 
by  the  rate  at  which  fresh  supplies  of  oxygen  can  be  transported  to  the  char  surface. 

When  a  laminar  boundary  layer  exists  over  the  char  surface,  the  transport  mechanism  of 
interest  is  molecular  diffusion  and  this  maximum  rate  is  referred  to  as  the  diffusion  - 
controlled  oxidation  rate.  Since  carbonaceous  materials  possess  extremely  high  intrin¬ 
sic  chemical  reactivities  at  temperatures  above  3000°  R  (1667°  K)  (refs.  1  and  2),  the 
char -layer  oxidation  is  generally  assumed  to  be  diffusion  controlled  under  most  reentry 
conditions.  (See  refs.  1  to  5.) 

In  the  present  report,  an  analytical  study  of  the  diffusion-controlled  oxidation  of  an 
ablative  char  layer  is  presented  and,  within  a  range  of  pressure  and  temperature  which 
is  typical  of  stagnation-point  conditions  during  orbital  reentry,  an  approximate  closed- 
form  solution  is  obtained.  Although  this  solution  is  based  on  a  number  of  simplifying 
assumptions,  it  is  believed  to  be  a  sound  conceptual  picture  of  the  diffusion-controlled 
oxidation  process  which  demonstrates  the  role  of  pyrolysis  gases  in  determining  the 
oxidation  rate  (and  hence  the  steady -state  pyrolysis  rate)  and  provides  a  relation  between 
the  original  chemical  composition  and  the  cold-wall  effective  heat  of  ablation  of  a 
charring  ablator. 


SYMBOLS 

All  dimensional  quantities  are  presented  in  terms  of  both  U.S.  Customary  Units  and 
the  International  (SI)  System  of  Units.  (See  ref.  6.) 

A,  B,C,D  stoichiometric  coefficients  per  mole  of  boundary -layer  gas  mixture 

(see  eq.  (100)) 

a,b, .  .  .,o  equilibrium  coefficients  (see  eqs.  (71)  to  (85)) 

cp  specific  heat  of  boundary -layer  gas  mixture,  British  thermal  units  per 

pound-°R  (joules  per  kilogram-°K) 
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coefficient  of  diffusion,  feet2  per  second  (meters2  per  second) 

mass  of  char  produced  by  decomposition  of  unit  mass  of  virgin  plastic 

stream  enthalpy,  British  thermal  units  per  pound  (joules  per  kilogram) 

enthalpy  potential  across  boundary  layer,  He  -  Hw,  British  thermal  units 
per  pound  (joules  per  kilogram) 

overall  heat  of  reaction,  British  thermal  units  per  mole  of  carbon  (solid) 
(joules  per  mole  of  carbon  (solid))  (see  eqs.  (113)  and  (116)) 

heat  of  formation  of  species  i,  British  thermal  units  per  mole  of  i  (joules 
per  mole  of  i) 

convective  heat-transfer  coefficient,  pounds  per  foot^-second  (kilograms 
per  meter^ -second) 

mass-transfer  coefficient,  pounds  per  foot2 -second  (kilograms  per 
meter2 -second) 

mass  fraction  of  species  i,  mass  of  i  per  unit  mass  of  gas  mixture 

elemental  mass  fraction,  mass  of  chemical  element  i  (no  matter  what 
chemical  form  it  may  appear  in)  per  unit  mass  of  gas  mixture 

mass  of  species  i  per  unit  mass  of  gases  injected  into  boundary  layer 

equilibrium  constant  (expressed  in  terms  of  partial  pressures)  for 
reaction  presented  in  equation  (n) 

thermal  conductivity  of  boundary -layer  gas  mixture,  British  thermal  units 
per  foot-second-°F  (joules  per  meter-second-°K) 

exponent  in  transpiration  cooling  expression  (see  eqs.  (125)  and  (127)) 

molecular  weight  of  gas  mixture  at  char  surface 

average  molecular  weight  of  gases  injected  into  boundary  layer 


Mi  molecular  weight  of  species  i 

m  mass  rate  of  oxidation  or  ablation,  pounds  per  foot2 -second  (kilograms 

per  meter2-second) 


r Uq  mass  of  carbon  per  unit  mass  of  virgin  plastic 

mjj  mass  of  hydrogen  per  unit  mass  of  virgin  plastic 

nij^  mass  of  nitrogen  per  unit  mass  of  virgin  plastic 

niQ  mass  of  oxygen  per  unit  mass  of  virgin  plastic 


na 

NLe 

nik 

P 

Pi 

Q* 

Qeff 


q 


^CW 


AvogadroTs  number 

t  •  u  ^D12cP 

Lewis  number,  — L — — 

K 

number  of  moles  of  species  i  per  mole  of  species  k  (see  eq.  (23)) 

pressure,  pounds  per  foot2  or  atmospheres  (newtons  per  meter^) 

partial  pressure  of  species  i,  pounds  per  foot2  (newtons  per  meter2) 

thermochemical  heat  of  ablation,  British  thermal  units  per  pound  (joules 
per  kilogram)  (see  eq.  (142)) 

q  ° 

cold-wall  effective  heat  of  ablation,  ~C'— - ,  British  thermal  units  per 

^vp 

pound  (joules  per  kilogram) 

convective  aerodynamic  heat  flux,  British  thermal  units  per  foot2 -second 
(joules  per  meter2-second) 

cold-wall  aerodynamic  heat  flux,  British  thermal  units  per  foot2 -second 
(joules  per  meter2 -second) 


qnet  net  heat  flux  to  surface  of  ablator,  used  in  defining  Q*,  British  thermal 

units  per  foot2-second  (joules  per  meter2 -second)  (see  eq.  (142)) 
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heat  flux  due  to  both  convection  and  chemical  reactions,  British  thermal 
units  per  foot^-second  (joules  per  meter^- second) 

r  nose  radius  of  hemispherical  body,  feet  (meters) 

T  temperature,  °R  (°K) 

v  boundary-layer  velocity  normal  to  char  surface,  feet  per  second  (meters 

per  second) 

W  ratio  of  diffusive  oxygen  flux  to  char  oxidation  rate,  defined  by 

equation  (95) 

Xj[  mole  fraction  of  species  i,  moles  of  i  per  mole  of  gas  mixture 

Xi  moles  of  species  i  per  mole  of  pyrolysis  gas 

y  distance  normal  to  char  surface,  positive  when  directed  away  from  char, 

feet  (meters) 

rj  transpiration  cooling  factor  (see  eq.  (120)) 

A  ratio  of  virgin  plastic  pyrolysis  rate  to  char  oxidation  rate  (see  eq.  (7)) 

p  boundary -layer  density,  pounds  per  foot^  (kilograms  per  me  ter  3) 

<j,T,v,cp  indices  in  empirical  polymer  formulas,  indicating  the  respective  number 

of  carbon,  hydrogen,  nitrogen,  and  oxygen  atoms  in  the  plastic 


Subscripts: 

C,H,N,0  pertaining  to  carbon,  hydrogen,  nitrogen,  and  oxygen,  respectively 

c  pertaining  to  ablative  char  layer 

ch  pertaining  to  or  resulting  from  chemical  reactions 

e  evaluated  at  outer  edge  of  boundary  layer 
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ij,k  pertaining  to  ith,  jth,  or  kth  chemical  species  or  element 

s  evaluated  at  stagnation  conditions 

vp  pertaining  to  virgin  plastic 

w  evaluated  at  char  surface 

12  denoting  values  for  binary  gas  mixture 

Superscripts: 

o  evaluated  for  zero  mass  injection  into  boundary  layer 

*  steady -state  solution 


ANALYSIS 


In  this  section,  the  governing  equa¬ 
tions  and  the  assumptions  used  in  the  pres¬ 
ent  analysis  are  stated,  and  the  ranges  of 
temperature  and  pressure  within  which  an 
approximate  closed-form  solution  is 
obtainable  are  identified.  The  distance  y 
used  in  this  analysis  and  the  various  zones 


^  X 

Mass  fraction!  *  Boundary-layer 
boundary  I  y  species 
layer 

Char  i 


r — Outer  edge  of 
/  boundary  layer,  e 


Char  surface  (w) 


Pyrolysis  interface 


Virgin 

P'aS,icl  •  ■  L-64-10. 940.1 

Figure  L-  Charring  ablator  showing  zones  and  interfaces. 


and  interfaces  pertinent  to  the  behavior  of  a  charring  ablator  as  analyzed  herein  are  pre¬ 
sented  in  figure  1.  The  virgin  plastic  is  assumed  to  contain  only  nitrogen,  oxygen,  carbon, 
and  hydrogen. 


Governing  Equations 

Consider  the  decomposition  of  a  hydrocarbon  polymer  according  to  the  equation 

CaHg-N^O^  -  Pyrolysis  gas  +  C(s)  (1) 

At  the  char  surface,  mass  conservation  of  the  element  nitrogen  may  be  represented  by 
the  equation 
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,n(kn, 


e‘%, 


%pv) 


Diffusion  away  from 
char  surface 


Convection 
away  from 
char  surface 


;MC  +  rMH  +  +  cpMq 

Flow  of  pyrolysis  gases 
from  pyrolysis  zone  to 
char  surface 


where  a,  r,  v,  and  <p  represent  the  relative  number  of  atoms  of  carbon,  hydrogen, 
nitrogen,  and  oxygen  in  the  virgin  plastic,  and  hm  ^  is  the  elemental  mass -transfer 
coefficient  for  nitrogen,  which  is  defined  by  the  equation 


hm,NiKN,e  KN, 


The  quantity  in  parentheses  on  the  right-hand  side  of  equation  (2)  represents  the  mass  of 
elemental  nitrogen  contained  in  the  pyrolysis  gases  produced  by  decomposition  of  a  unit 
mass  of  virgin  plastic. 

In  order  for  equation  (3)  to  be  rigorously  valid,  all  the  species  containing  nitrogen 
must  have  nearly  equal  diffusion  coefficients.  Since  subsequent  equations  will  place 
similar  requirements  on  the  species  containing  oxygen,  hydrogen,  and  carbon,  the  use  of 
the  diffusion  coefficients  D^,  D q,  Djj,  and  D^.  implies  that  the  diffusion  coefficients 
of  all  the  species  are  nearly  equal  -  that  is,  the  gas  mixture  is  effectively  binary.  This 
binary  mixture  requirement  is  discussed  in  the  section  "Review  of  Assumptions."  Equa¬ 
tions  similar  to  equation  (2)  may  be  written  for  oxygen,  carbon,  and  hydrogen  as  follows: 


-hm  o(^0  e  '  w)  +  =  1 - - pvp 

m,OV  H},e  H),w;  V  ^  h  Lm  +  tMh  +  +  cpMQ  Vp 


■hm,c(KC,e  -  Kc,w)  +  (&C Pv)w  =  (— - — - - - —  *  ™vp  +  mc  (5) 

x  7  v  7  \aMc +  tMh  +  +  ^Mo  / 


c  +  tMjj  +  i/Mn  +  cpM.Q 


where  f  represents  the  mass  of  char  produced  by  decomposition  of  a  unit  mass  of 
virgin  plastic.  Now  define  X  as 
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and  note  that 


(pv)w  =  mc  [l  +  A(1  -  f)_ 


(8) 


The  respective  masses  of  nitrogen,  oxygen,  carbon,  and  hydrogen  per  unit  mass  of  virgin 
plastic  can  be  expressed  as 


I'M 


mN  “ 


m0  = 


m^  = 


mu  = 


N 


H 


oMq  +  tMh  +  i/Mj^  +  <pM0 

_ _ 

oMc  +  tMh  +  i/Mn  +  (j)Mq 

_ _ 

oMq  +  tMh  +  nMN  +  0Mq 

_ _ 

+  tMh  +  I'Mjq-  +  i^Mq 


(9) 


(10) 


(ID 


(12) 


Substitution  of  equations  (7)  to  (12)  into  equations  (2),  (4),  (5),  and  (6)  yields  the  following 
expressions  for  the  elemental  mass  fractions  at  the  char  surface  /  by  definition, 

T 7"  T  7-  r\  \  \ 


KC,e  =  K 


H,e  -  °): 


Z-  KN,ehm,N  +  mNAlilc 

%,w  - - 2 - 2 - 


h 


m.N  +  C1  +  A^  -  f3 


(13) 


*0, 


W 


KQ,ehm,Q  +  mOAlilc 

hm,0  +  L1  +  -  f3 


(14) 


K 


C,w 


1  +  (mc  -  f] 

£] 

mc 

hm, 

C  + 

[l  +  X(1 

0] 

rii 

(15) 


K 


mjjXmc 


H,w 


hm,H  +  C1  +  A^  -  f3™< 


(16) 
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Equations  (13)  to  (16)  are  the  governing  equations  for  this  study.  These  equations, 
however,  cannot  be  solved  unless  the  mass-transfer  coefficients  hm^  can  be  evaluated. 


Evaluation  of  Mass -Transfer  Coefficient 

In  this  report,  the  mass -transfer  coefficients  hm  ±  are  expressed  in  terms  of 
the  convective  heat-transfer  coefficient  through  the  application  of  an  analogy  between  the 
transfer  of  energy  and  mass  in  a  laminar  boundary  layer.  This  analogy  (with  minor 
variations)  is  often  used  in  treating  the  diffusion -controlled  oxidation  of  graphite.  (See 
refs.  2,  3,  7,  8,  and  9.) 

Lees  in  reference  9  shows  that,  for  the  case  of  a  laminar,  unit  Lewis  number, 
stagnation -point  boundary  layer,  the  ith  elemental  mass  fraction  and  the  total  enthalpy 
are  related  by  the  following  equation: 


Ki  -  Ki,w 


HS  ^s,w 


Hence 


Ki.e  "  Ki.w  Hs.e  ~  H: 


Kj,e  -  Ki,w  9Hs 

9y  "Hs.e-Hs.way 


and  since 


PD12c, 


it  follows  that 


Pd12 


Ki,e  -  Ki)W  NLe 


^s.e  -  ^s,w 


^i,e  “  ^i/w)*1 


But  the  ith  elemental  mass -transfer  coefficient  is  defined  by  the  equation 


pDi  fa~=  (Ki’e  ’  Ki,w)hm,i 


Hence,  if  the  boundary -layer  mixture  is  effectively  binary, 


^m,i  “  ^ 


This  assumption  is  discussed  in  the  section  "Review  of  Assumptions." 
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Pyrolysis  Gas  Composition,  Char  Composition,  Assumed  Reactions, 
and  Equilibrium  Expressions 

The  elemental  mass  fractions  which  appear  in  equations  (13)  to  (16)  can  be 
expressed  in  terms  of  the  mole  fractions  of  the  various  chemical  species  by  the  following 
equation: 


where 


Mi  atomic  weight  of  elemental  species  i 
M  molecular  weight  of  gas  mixture 

Xk  mole  fraction  of  kth  chemical  species  in  gas  mixture 
nik  number  of  moles  of  i  per  mole  of  k 


(23) 


The  chemical  species  which  must  be  considered  are  determined  by  the  reactions  which 
take  place  between  the  pyrolysis  gases,  the  boundary -layer  species,  and  the  char. 

The  elemental  mass  fractions  Kq  w  are  expressed  in  terms  of  the  reactants  and 
products  under  the  assumption  that  chemical  equilibrium  (both  heterogeneous  and  gas 
phase)  exists  at  the  char  surface.  The  pyrolysis  gases  are  assumed  to  be  composed  of 
Ng,  H 2>  NHg,  CO,  CO2,  H^O,  ^2^2’  ^2^4’  anc*  anc*  c^ar  *s  assumed  to  be 
graphite. 

The  following  chemical  reactions  are  considered: 


2NH3  z  3H2  +  N2 

(24) 

CH4  +  202  z  C02  +  2H20 

(25) 

C2H2  +  |  °2  -  2C02  +  H2° 

(26) 

C2H4  +  202  -  2CO  +  2H20 

(27) 
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l°2=° 

(28) 

|H2  -H 

(29) 

|h2  +  1o2=oh 

(30) 

H2  +  J  °2  :  H2° 

(31) 

|N2  -N 

(32) 

iN2  +  i°2rN° 

(33) 

C(s)  Z  C(g) 

(34) 

C(s)  +  io2  Z  CO 

(35) 

C(s)  +  C02  z  2 co 

(36) 

c(s)  +  h2o  ;h2  +  co 

(37) 

C(s)  +  N  z  CN 

(38) 

C(s)  +  NH3  z  HCN  +  H2 

(39) 

The  notation  s  or  g  indicates  species  in  a  solid  or  gaseous  state,  respectively.  The 
equilibrium  constants  for  reactions  (eqs.  (24)  to  (39))  may  be  written  as  follows: 


PH23pN2  S3XN2  2 

HP, (24)  n  2  X  2  P 


'NH, 


NH, 


PC0oPHo0  XC0oXHo0 


K 


2  2 


p’<25) "  Pch4Po22 


xch4xo2 


(40) 


(41) 
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(42) 


K 


PC022pH2Q  XC022xH2Q  _1/2 


P’(26)  PC  H  p0  5/2  XC  H  X0  5/2 
u2  '•'2il2  u2 


D  2n  2  X  2X  2 

pco  ph9o  aco  aH0O 
K  -  —  = - p 


^P,  (27) 


pC2H4P02^  XC2H4X02 


(43) 


y  _  P0  _  X0  „l/2 

W^‘V^P 


(44) 


K  =  PH  XH  1/2 

P’^29^  pH  1/2  xH  V2 
H2  H2 


(45) 


K 


P0H  _  X0H 
P>(30)  _  1/2  1/2  "x  l/2Xn  1/2 

n2  u2  H2  U2 


(46) 


K 


0h2O  XH2°  -1/2 


P'<31)  PH  p0  1,2  XH  XQ  1/2 
n2  u2  h2  U2 


(47) 


Y  _  PN  _  XN  1/2 

p’(32r^'V^p 


(48) 


K 


pNO 


NO 


P,(33)  -  ^^1/2^172  -  XnV2XqV2 


(49) 
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Kp,(34)-pC(g)-XC(g)p 


(50) 


_  PCO  _  XCO  1/2 
P,(35)  '  ^  1/2  -  x^l/2 


Pco2  xco2 

?,(36)  -  pc02  -  xco 


ph2pco  xh2xco 

Kp,(37)-  Ph2o  "  p 


PCN  _  XCN 
^p,(38)  pN  xN 


PHCNPH2  xhcnxh2 


p,  (39)  P;NH, 


Equations  (40)  to  (55)  may  now  be  used  to  express  the  mole  fraction  of  each  species  in 
terms  of  the  mole  fractions  of  Ng,  Hg,  and  CO.  The  resulting  expressions  are 


X  -ax  3/2X  1/2 

xnh3  _  ^Hg  an2 


XCH4  "  bXH2' 


Xp  TT  "  CXtt 
^2^2  W2 


13 


xr  R  =  dX  ' 

-^2 


xo  -  eXco 


XH  =  «H. 


1/2 


XOH  _  gXCOXH, 


1/2 


xH2o  =  ^h/co 


XN  =  «N. 


1/2 


X  —  iX 

ano  _  ]AN2  aco 


XC(g)  =  k 


xo2  _  zxco” 


xco2  _  mXco 


XCN  “  nXN21//2 


X  -  oX  ^/^x 

ahcn  -  0AH2  an2 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 
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where 


P 


KP,(24)1/2 


(71) 


b  = 


KP,(35)  P 


Kp,(37)2Kp,(36)Kp,(25) 


(72) 


K 


c  = 


Pj  (35) 


KP>  (36)  KP,(37)Kp,  (26) 


(73) 


d  = 


KP,(35)  p 


KP>(27)Kp,  (37) 


(74) 


e  = 


^,(28) 


K 


P,  (35) 


(75) 


f  =  K 


p,(29)p 


-1/2 


(76) 


g  = 


KP,(30)p 


1/2 


K 


P,(35) 


(77) 


h  = 


P 


K 


1  =  KP,(32)P 


P,(37) 

-1/2 


(78) 


(79) 
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(80) 


J  = 


KP,(33)p 


1/2 


K 


p,(35) 


k  =  K 


p,(34)p 


-1 


(81) 


l  = 


K  2 
Pj  (35) 


(32) 


m  = 


P 


K 


P,(36) 


(83) 


n  —  K .  /nn\K 


P>  (38)  p,(32)p 


-1/2 


(84) 


K 


o  = 


P,(39) 


KP,(24)1/2 


(85) 


The  governing  equations  (eqs.  (13)  to  (16))  may  now  be  restated  as  follows: 


M 


- _ 5/gY  3/2«-  1/2  1/2^-  ,  1/2  -Y  1/2 

N,w  ^  pH2  XN2  +  ^XN2  XCO  +  nXN2  +  lXN2  +  2XN2 


♦  0X„  V2X  l/2\  .  *  mNXli,c 


(86) 


=  1/2- 


+  “h^co  +  x 


CO 


w 


V  +  m0Amc 
h  +  [l  +  A(1  -  fj]  mc 


(87) 
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XC,W  -  ^K„2  ♦  XC0  +  2cXH„  *  k  +  “^CO2  +  „1/2  +  °XH„1/2xN.  1/2 


M  V  2 


2  2 


+  2dX 


H. 


w 


e- 

(mc  - f) 

“c 

h  + 

[l  +  X(1  -  fj 

]“c 

(88) 


XH,w  '  -r(3aXH23/2xN21/2  +  4bXH22  +  2XH2  +  2cXH2  +  “H,172  +  6XCOXH21/2 


mHXmc 


+  2hX„  Xrn  +  oXtt  1/2Xn  1/2  +  4dX„  2  _  _ 

H2  CO  H2  N2  H2  jw  h  +  [l  +  X(1  -  f )] 


(89) 


nr 


Approximate  Solution 

General  expressions.-  When  all  the  equilibrium  coefficients  (a,  b,  etc.)  are  «1.0, 
equations  (86)  to  (89)  may  be  approximated  as  follows: 


K.t  ...  =  2  ...  -  K 


KN,eh  +  mNXmc 


N’w  M  N2’w  N2,w  h  +  [l  +  X(1  -  f )]  m 


(90) 


~  M0 

K°’w  =  ^Xc°’w 


M 


O 


K, 


K0,eh  +  mOXlilc 


MCO  CO’W  h  +  C1  +  x(!  "  f)]m( 


(91) 


M. 


M, 


K 


C,w-  jj  CO,w 


K. 


M 


CO 


[1  + 

(mc  -  f)x] 

“c 

h  + 

[l  +  X(1  -  f) 

]”c 

(92) 


2M 


K 


H 


H,w 


Xtt  =  K 


mTTXni 
H  c 


M 


H2’w  H2’w  h+  [l  +  X(1  -  f )]  m( 


(93) 
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and  equations  (91)  and  (92)  may  be  solved  for  mc  as  follows: 


where 


(94) 


W  = 


-  mOX 


(95) 


It  should  be  noted  that  the  solution  for  mc,  presented  in  equations  (94)  and  (95),  is  valid 
for  the  general  case  of  transient  ablation.  The  parameter  W  is,  however,  a  function 
of  A  and  in  the  transient  case  A  is  not  known  but  must  be  obtained  from  experiment 
or  through  the  iterative  solution  of  a  transient  ablation  computer  program.  For  the 
steady-state  ablation  case,  however,  A  is  a  known  constant.  The  steady-state  ablation 
case  is  considered  in  detail  in  the  subsequent  section  of  this  paper. 

The  boundary -layer  composition  at  the  char  surface  can  be  obtained  from  equa¬ 
tions  (90),  (91),  and  (93)  as  follows: 


K, 


W^+mNX 

Ko.e  N 


K 


CO,w 


K 


H2’w  W 


w  + 

1  +  A(1  - 

f) 

K0,e 

MCO 

[w  +  h1q> 

M0 

] 

w  + 

1  +  A(1  - 

■  f) 

K0,e 

mHA 

K 


+  1  +  A(1  -  f) 


0,e 


(96) 


(97) 


(98) 


The  mass  fractions  K,  can,  of  course,  be  expressed  in  terms  of  mole  frac- 

Ij  w 

tions  X-  through  the  relation 
Ij  w 
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(99) 


5. 

Mi 

Xi  = - 1 


yK- 

Z  M 


Once  the  boundary -layer  composition  at  the  char  surface  has  been  determined,  an 
overall  heat  of  reaction  may  be  calculated  as  follows.  An  equation  for  the  overall  reac¬ 
tion  per  mole  of  boundary-layer  gas  mixture  may  be  written  as 


[A02  +  BN2 


X 


CO, 


C  CO  + 


co„  + 


CO 


2  X 


Xprr  X^i  rr  X^i  IT 

CXI,  Onlln  Onll  , 

•  CH.  + - ^  CoHr,  +  —Z  -  CnE 


CO 


Y  2  2  Y 

ACO  ACO 


2  4 


‘N9  XNH9  XH0O  xh9 

—  N9  +  — — -  NH9  +  H90  +  — — —  H9 

“  XT'  O  xr  Ci  Y7"  ^ 


CO 


CO 


CO 


CO 


[DC(s)] 


XN2,wN2  +  XCO,wCO  +  X 


H2,wH2j 


(100) 


where  X^  is  the  number  of  moles  of  species  i  per  mole  of  pyrolysis  gas.  The  terms 
of  equation  (100)  are  identified  as  follows: 


^Species  from  stream]  +  [Pyrolysis  gasesj  +  jchar]  Z.  [Reaction  products] 


The  requirement  that  the  number  of  atoms  of  a  given  element  must  be  the  same  on  both 
sides  of  equation  (100)  results  in  the  following  relations: 


2A  + 


C  1  + 


2xco2  %2o' 

XCO  XCO, 


=  X 


CO.w 


/ 


2B 


2XN 


2’ 


w 


(101) 


(102) 
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xco„  xch4  2xc2h2 

C  1 1  +  — — —  + - -  + - — + 


2X 


C2H4 


+  D  =  X 


xco  xco 


X 


CO 


X 


CO 


CO,w 


(103) 


,4XCH4  2XC?H9  4XC?H4  3XNH,  2XHpO  2XHp 
Cl  — - 4+  2  2  +  -_A  -1  +  — - -  +  — r— ^ —  +  — — - 1  =  2Xtt  w  (104) 

x~  xco  xco 1  2’ 


CO 


X 


CO 


CO 


CO 


2XC02  xh2o 

but,  for  each  mole  of  reaction  products  produced,  NAC  II  +  — - -  +  — - )  oxygen 

A  i  v  y 

ACO  ^CO 

atoms  entered  the  reaction  in  the  form  of  pyrolysis  gases  and  N^D  atoms  of  carbon 
entered  the  reaction  in  the  form  of  char.  Hence,  from  the  definition  of  f  and  niQ, 


2XCO?  XHpO 
1  +  — - -+-=-£ 


\  Xco  Xco  /  _  m0McX 


D 


M 


O 


(105) 


Similarly, 


X 


CO, 


C  1  + 


+  . 


X 


CO 


D 


=  (mc  -  f)x 


(106) 


Ax 


CH 


+  . 


X 


CO 


D 


mH^C^ 

MH 


(107) 


2X 


N, 


+  . 


CO 


D 


\ 


/  _  mNMCX 


M 


N 


(108) 
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and  equations  (101),  (102),  (103),  and  (104)  with  the  aid  of  equations  (105)  to  (108)  can  be 
expressed  per  mole  of  C(s)  (char)  by  dividing  through  by  D.  Hence, 


2A  mOMCx  XCO,w 
D  M0  “  D 


(109) 


2B  ,  mNMCX 
D  Mn 


2Xn2’W 

D 


(110) 


V  -  f)x  +  l  = 

'  u  /  D 


(HD 


mH^C^  ^2,w 

mH  D 


(112) 


Since  the  overall  heat  of  reaction  is  defined  as 

ahr  =  Enthalpy  of  products  -  Enthalpy  of  reactants 

from  equation  (100),  the  overall  heat  of  reaction  per  mole  of  solid  carbon  (char)  consumed 
may  be  stated  as 


AHR  Cq,W  AHf , CO  '  jJi  (XCOAHf, CO  +  XC02 AHf , C02  +  XCH4AHf , CHj 


+  XC2H2AHf,C2H2  +  XC2H4AHf,C2H4  +  XNH3AHf,NH3  +  XH20AHf,H20 


(113) 


but  according  to  equations  (111)  and  (105) 


21 


(114) 


XCO,w 

D 


(mC  -  f)x  +  1 


C 

DXco 


m0McX _ 1 _ 

Mo  xco  +  2XC02  +  xh2o 


(115) 


Hence  equation  (113)  becomes 


AH 


R 


(mc  -  f)x  +  1  AH^cq 


mOMC^ 


I5* 


aHf.i 


M 


0  xco  +  2xco2  +  xh2o 


(116) 


The  heat  flux  due  to  chemical  reactions  (recombination  reactions  for  oxygen  and  nitrogen 
are  not  included  since  they  are  usually  accounted  for  in  computing  convective  heat  fluxes) 
may  now  be  expressed  as 


^ch  “  “ 


AHRmc 

MC 


hK 


0,e 


MCW 


4mc  -  f)x  +  i]AHf,co  " 


mQM^X 

Mq 


I*. 


AH 


xco  +  2xco2  +  xh20 


(117) 


(The  minus  sign  is  used  in  equation  (117)  since,  by  definition,  an  exothermic  reaction  has 
a  negative  heat  of  reaction.)  Equation  (117)  may  now  be  used  to  define  an  "effective 
chemical  enthalpy  potential  AHch"  as  follows: 


AHch  .  M 

cn  h  MpW 


I  *iAHf,i 

Jr,  >  1  m0MrX  4 

f"c  ■  ^ +  1]4Hf"co '  IT73 


I 


CO  T  “"xco2  +  xh2o 


>  (118) 


Hence,  the  total  heat  flux  (convective  +  chemical)  is  expressible  as 


qt  =  h(AH+AHch) 


(119) 
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Since  expressions  for  the  gas  composition  at  the  char  surface  have  been  derived, 
consideration  of  the  transpiration  cooling  effect  produced  by  the  pyrolysis  gases  is  of 
interest.  In  practical  applications  the  reduction  in  heat  transfer  which  results  from  the 
injection  of  gas  into  the  boundary  layer  is  often  accounted  for  by  using  the  relation 

h  =  h°  -  17 (pv)w  (120) 


and  in  reference  10  it  is  suggested  that  the  transpiration  factor  17  varies  approximately 
as  M"^  where  M  is  the  molecular  weight  of  the  injected  gas  and  l  ranges  from 
1/3  to  1/4.  If  the  gas  phase  combustion  reactions  all  occur  near  the  char  surface,  the 
species  which  are  actually  injected  into  the  boundary  layer  are  the  combustion  products, 
that  is,  Ng,  Hg,  and  CO.  The  mass  fraction  of  Ng  in  the  injected  gases  is  expressible  as 


and  similarly 


m 


N, 


K 


N, 


mNmvp 


mNX 


(pv) 


w 


(1  -  f)mvp  +  mc  (1  -  f)A  +  1 


(121) 


m 


H, 


K 


H, 


w 


mTTm 
H  vp 


mHX 


(1  -  f)mvp  +  mc  (1  -  f)A  +  1 


(122) 


K 


CO 


_  “CO  (mO  +  mC  -  f)“vp  +  “c  _  (mO  +  mc  -  f)A  +  1 


(pv) 


w 


ZE 

(!  -  f)myp  +  m( 


(1  -  f)A  +  1 


(123) 


Hence,  the  average  molecular  weight  of  the  injected  gases  is  given  by  the  expression 


M  = 


(1  -  f)A  +  1 


mNX  ,  mHX  , 


+  m^-,  -  f  jA  +  1 


M 


N, 


M 


H, 


M 


CO 


(124) 


Now  the  ratio  of  the  transpiration  cooling  factor  for  a  gas  of  arbitrary  molecular  weight 
to  the  transpiration  cooling  factor  for  carbon  monoxide  (M  =  28)  is 
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^M=28 


/mV* 

(1  -  f)X  +  1 

-i 

lSm^A 

1  i  ** 

\  28  / 

+  14m^X  +  itIqA  +  m^A  -  f A  +  1 

(l-f)X+l 

(125) 


and  since  for  a  laminar  boundary  layer 

r?M=28  ~  ^air  ~  0,7 

equation  (120)  may  be  approximated  by  the  following  expression 

l 


(126) 


h  «  h°  -  0.7 


1  + 


13mjjX 


(1  -  f)A  +  1 


(pv) 


w 


(127) 


Steady-state  expressions.-  The  solutions  obtained  in  this  analysis  (as  presented  in 
eqs.  (94)  to  (98),  (116)  to  (119),  and  (127))  are  applicable  to  the  case  of  transient  ablation  - 
that  is,  the  case  where  the  char  surface  and  the  pyrolysis  interface  move  at  different 
linear  rates.  However,  these  solutions  depend  on  the  ratio  of  pyrolysis  rate  to  char 
oxidation  rate  X,  and  X  is  not  known  a  priori  except  in  the  case  of  steady-state  abla¬ 
tion  where  X  =  1/f.  The  limiting  case  of  steady  ablation  is,  however,  an  important  one 
and  hence  the  solutions  for  X  =  1/f  are  presented  in  the  following  equations: 


„  Kn  h 

—  =  m  =  -Qig- 
vp 


(128) 


where 


Mn 

fW*  = - mr  -  = 

M„  ^  u 


M0(ct  -  <p) 


oMq  +  tMh  +  fMN  + 


(129) 


K 


N2,w 


(130) 
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(131) 


h  -  h°  -  r?mvp  =  h°  -  0.7  (l  +  UmRj  myp  (135) 

DISCUSSION 

In  this  section,  the  assumptions  used  in  the  analysis  are  reviewed,  a  physical 
interpretation  of  the  approximate  solution  is  presented,  and  the  effects  of  virgin  plastic 
composition  and  pyrolysis  gas  composition  on  the  char  oxidation  rate,  the  steady-state 
ablation  rate,  the  cold-wall  effective  heat  of  ablation  (defined  subsequently),  and  the  heat 
release  due  to  chemical  reactions  are  discussed.  In  the  section  "Analysis/'  both  general 
and  steady-state  expressions  of  the  approximate  solution  were  presented.  Before  a 
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solution  to  the  general  expressions  can  be  obtained,  the  ratio  of  pyrolysis  rate  to  char 
oxidation  rate  X  (which  is  not  known  except  for  the  steady-state  ablation  case)  must  be 
specified.  Hence,  while  the  general  expressions  can  be  used  in  conjunction  with  a  com¬ 
puter  program  (such  as  that  presented  in  ref.  11)  to  yield  transient  ablation  rates,  they 
do  not  lend  themselves  to  demonstrating  the  general  characteristics  of  diffusion- 
controlled  char  oxidation.  For  steady -state  ablation,  on  the  other  hand,  X  =  1/f  and  as 
a  result  the  equations  are  simple  and  provide  a  clear  physical  picture  of  diffusion- 
controlled  char  oxidation.  Accordingly,  the  following  discussion  is  restricted  to  the 
steady-state  case.  Any  remarks  concerning  the  char  oxidation  rate  apply  equally  well 
to  the  virgin -plastic  ablation  rate  since  mc  =  fri^  . 

Review  of  Assumptions 

The  major  assumptions  used  in  the  present  analysis  are: 

(a)  The  virgin  material  is  composed  of  only  oxygen,  hydrogen,  nitrogen,  and 
carbon. 

(b)  The  mass-transfer  coefficients  for  all  gas  phase  species  are  numerically  equal 
to  the  convective  heat -transfer  coefficient. 

(c)  All  the  equilibrium  constants  in  equations  (86)  to  (89)  are  «1.0. 

(d)  Chemical  equilibrium  (both  gas  phase  and  heterogeneous)  is  obtained  at  the  char 
surface. 

(e)  The  transpiration  factor  77  varies  inversely  with  the  kth  power  of  the  mean 
molecular  weight  of  the  gases  injected  into  the  boundary  layer. 

(f)  The  pyrolysis  gases  are  composed  of  N2,  Hg,  NHg,  CO,  C02,  C^H^,  C2H2, 

CH^,  and  H20. 

Most  presently  used  high -temperature  resins  (silicone  resins  and  rubber  are 
notable  exceptions)  are  composed  of  O,  H,  and  C.  Nylon,  which  is  commonly  used  as  a 
reinforcement  in  ablative  materials,  is  composed  of  O,  H,  N,  and  C;  hence,  mixtures 
such  as  epoxy  nylon  or  phenolic  nylon  can  be  treated  by  the  techniques  presented  herein. 
Also,  the  present  method  could  be  applied  to  a  material  composed  of  a  hydrocarbon  resin 
and  an  inert  filler. 

The  assumption  made  in  evaluating  the  mass -transfer  coefficients  is  commonly 
made  in  analyses  of  carbon  combustion.  (See  refs.  2,  3,  7,  8,  and  9.)  In  reference  9 
this  assumption  is  shown  to  be  valid  whenever  the  gas  mixture  can  be  considered  binary 
with  respect  to  diffusion  and  the  Lewis  number  is  unity.  Even  though  neither  of  these 
requirements  are  rigorously  met  for  carbon  combustion,  use  of  the  mass-transfer- 
coefficient  assumption  yields  results  which  are  in  qualitative  agreement  with  experiment. 
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A  measure  of  the  validity  of  this  assumption  for  carbon  combustion  can  be  obtained 
by  comparing  the  present  solution  with  that  presented  in  references  1  and  4.  The 
analysis  presented  in  references  1  and  4  considered  multicomponent  diffusion  and  placed 
no  restriction  on  Lewis  number.  From  that  analysis  the  following  expression  for  the 
diffusion-controlled  stagnation -point  carbon  oxidation  rate  in  air  was  obtained  and  pre¬ 
sented  in  reference  4: 

=  6.35  x  10"3  lb/ft3/2-sec-atm1//2  =  0.05378  x  10'3  kg/m^-N^-sec  (136) 

Before  the  results  of  the  present  analysis  can  be  compared  with  equation  (136)  a  value 
for  the  heat -transfer  coefficient  must  be  obtained.  Evaluation  of  the  Fay  and  Riddell 
relation  (ref.  12)  has  shown  that  values  on  the  order  of 

3.92  x  10'2  lb/ft3/2-sec-atm1//2  =  0.033198  x  10"2  kg/m1/2-]*1/2 -sec  (137) 

are  typical  of  stagnation -point  conditions  on  a  manned  vehicle  during  entry  from  earth 
orbit.  If  this  value  and  the  oxygen  mass  fraction  for  air  (0.232)  are  substituted  into 
equations  (94),  (95),  and  (127)  there  results, 

=  6.07  x  10'3  lb/ft3/2-sec-atm1/2  =  0.051407  X  10'3  kg/m1/2-N1/2-sec  (138) 

which  agrees  with  equation  (136)  to  within  5  percent.  The  total  heat  flux  calculated  from 

equation  (119)  was  also  found  to  agree  with  the  values  presented  in  references  1  and  4 

to  within  5  percent.  Hence  the  unity  Lewis  number  assumption  may  be  considered  valid 

for  carbon  combustion.  Only  the  presence  of  light  gases  (primarily  Hg)  makes  the 

assumption  less  applicable  to  the  present  analysis  than  to  carbon  combustion.  Hence  the 

assumption  should  be  quite  good  for  low  values  of  free -stream  oxygen  mass  fraction 

( as  -  0  and  JL  -  0,  eq.  (98) )  and  although  the  validity  of  the  assumption 

\  u,e  w2?w  / 

decreases  with  increasing  Kq  ,  it  should  provide  an  acceptable  qualitative  description 
of  the  transport  process  for  air.  This  assumption  is  considered  in  the  light  of  calculated 
char-surface  boundary -layer  compositions  subsequently. 

The  approximate  solution  was  obtained  by  assuming  that  the  equilibrium  coefficients 
appearing  in  equations  (86)  to  (89)  were  small  compared  with  unity.  A  more  rigorous 
statement  of  this  requirement  would  be  that  the  sum  of  all  the  equilibrium  coefficients 
appearing  in  a  given  equation  must  be  small  compared  with  unity.  In  order  to  assess  the 
validity  of  this  assumption,  the  equilibrium  coefficients  were  evaluated  by  using  the  data 
presented  in  reference  16.  The  variation  of  these  coefficients  (evaluated  for  a  pressure 
of  1  atm)  with  temperature  is  presented  in  figure  2.  By  using  the  curves  presented  in 
figure  2,  ranges  of  char  surface  temperature  and  pressure  for  which  the  sum  of  the  equi¬ 
librium  coefficients  appearing  in  any  one  of  equations  (86)  to  (89)  was  <0.1  were  identified 
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The  ranges  of  char  surface  tem¬ 
perature  and  pressure  within  which 
the  sum  of  these  coefficients  is 


-f(H) 
k  (c (q  ) ) 
:c(c2h2) 

■O(HCN) 


d  (c2  H4: 

-b(CH.) 


<0.1  (and,  hence,  within  which  the 

n  (CN) 

approximate  solution  is  valid)  are 

presented  in  figure  3.  The  char  —  XwP'rX 

Jr  /^£clc2H2) 

surface  temperatures  and  pres-  o _ _ ^C^7^°(hcn) 

sures  experienced  by  a  manned  _  / 

ballistic  vehicle  (such  as  the  j  1  /  / /f  / 

Mercury  spacecraft)  during  entry  j  \1  /  / //  / 

from  earth  orbit  fall  within  this  P  \  /  / / /  /  / l(N) 

range  and  hence  the  present  _  "2p  \\  /  /  /  /  / 

approximate  solution  should  be  P  \v/  j  / 

applicable  to  such  vehicles.  On  o  -3j—  Vy\  <  / 

the  other  hand,  ballistic  missiles  I  L_  \\\y^  /  / 

and  lunar  return  vehicles  experi-  f  4  V  y(/\  /  / 

ence  temperatures  higher  than  ^  / \<  \  /\  /SC 

o  X  //\V  \  /-d<c2H4) 

those  shown  in  the  cross-hatched  dr  |  //  Nfv 

region  of  figure  3.  Solution  of  the  ‘5p"  y 

governing  equations  at  these  ~  h  //  yJ 

higher  temperatures  is  straight-  -6^—  1/  1  f  /X<h(H20} 

forward  but  would  require  the  use  —  /  i  I  P//  j(N0) 

of  a  digital  computer.  _7_  /  //  ^^m(co2) 

The  assumption  of  chemical  —  / 

equilibrium  at  the  char  surface  is  .6i _ i  .  I -i  .  I  _ J l 1  111  I  1 

.  -  0  1000  2000  3000  4000  5000  6000  °K 

a  natural  consequence  of  con-  ( _ | _ ( _ ( _ ^ _ 

sidering  the  case  of  diffusion-  0  1800  3600  5400  7200  9000 OR 

Temperature  ,  T 

Controlled  char  oxidation.  In  any  Figure  2.-  Variation  with  temperature  of  equilibrium  coefficients, 

kinetic  system,  the  slowest  step  Pressure,  l  atm  «U0l3MN/m2). 

is  rate  controlling.  Hence,  in  the  diffusion-controlled  case,  the  rate  of  transport  of 
reactants  must  be  slow  compared  to  the  rates  of  both  homogeneous  and  heterogeneous 
reactions.  Since  the  time  required  for  the  slowest  reaction  to  equilibrate  is  small  com¬ 
pared  with  the  time  required  to  deliver  fresh  quantities  of  reactants  to  the  char  surface, 
chemical  equilibrium  is  obtained. 


Temperature  ,  T 

Figure  2.-  Variation  with  temperature  of  equilibrium  coefficients. 
Pressure,  1  atm  (0.1013  MINI/ m2). 


The  assumption  made  with  regard  to  the  transpiration  factor  77  is  semiempirical 
in  nature  and  is  based  on  experiments  and  analyses  involving  nonreacting  boundary 
layers.  This  assumption  would  be  applicable  to  the  present  problem  if  all  the  gas  phase 
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reactions  occurred  in  a 
thin  zone  near  the  char 
surface.  When  this  sit¬ 
uation  occurs,  the  gases 
injected  into  the  bound¬ 
ary  layer  are  the  reac¬ 
tion  products  which 
diffuse  from  the  reac¬ 
tion  zone  through  the 
boundary  layer  without 
undergoing  further  reac¬ 
tion.  It  should  be  noted 
that,  in  analyzing  reentry 
heating,  several  investi¬ 
gators  have  assumed  a 
frozen  boundary  layer 
and  local  chemical  equi¬ 
librium  at  the  vehicle 
surface.  (See,  for 
instance,  ref.  1.)  This 
assumption  is  equivalent 
to  assuming  a  thin  reac¬ 
tion  zone  near  the  char 

surface.  When  reactions  occur  throughout  the  boundary  layer,  the  present  approach  is 
not  applicable;  however,  to  the  author Ts  knowledge,  no  solutions  for  the  transpiration 
effect  of  pyrolysis  gases  injected  into  a  reacting  boundary  layer  are  available.  Cer¬ 
tainly,  there  exists  the  possibility  that  transpiration  cooling  effects  significantly  different 
from  those  predicted  by  equations  (120)  to  (127)  will  be  obtained.  As  the  following  dis¬ 
cussion  develops,  however,  it  will  be  seen  that  transpiration  factors  which  differed  from 
the  presently  used  values  by  factors  of  two  or  even  three  would  not  significantly  change 
the  qualitative  results  of  this  analysis.  Hence  the  present  assumption  is  reasonably 
accurate  for  the  case  of  frozen  boundary  layers  with  chemical  equilibrium  at  the  char 
surface  and  could  in  no  case  introduce  errors  which  would  significantly  alter  the  con¬ 
ceptual  picture  of  diffusion-controlled  char  oxidation  presented  herein. 

The  pyrolysis  gases  assumed  in  the  present  treatment  represent  the  major  species 
which  have  most  often  been  reported  by  chemists  who  have  attempted  experimentally  to 
determine  pyrolysis  gas  compositions.  (See,  for  instance,  refs.  13  and  14.)  The  pyroly¬ 
sis  gases  considered  in  the  present  report  are  those  which  reach  the  char  surface. 

Hence,  while  the  initial  decomposition  of  the  plastic  may  produce  large  molecular  weight 
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hydrocarbon  species,  it  is  assumed  that  cracking  reactions  in  the  char  reduce  these 
large  fragments  to  the  species  listed  in  assumption  (f). 

Pyrolysis  gas  compositions,  calculated  for  the  equilibrium  decomposition  (i.e.,  all 
cracking  reactions  in  chemical  equilibrium)  of  phenolic  nylon  are  presented  in  refer¬ 
ence  15  and  contain  several  species  not  included  in  assumption  (f).  The  inclusion  of  all 
the  species  presented  in  reference  15  in  the  present  analysis  might  possibly  result  in 
small  changes  in  the  region  of  validity  shown  in  figure  3.  The  species  included  in 
assumption  (f)  do,  however,  account  for  approximately  90  percent  (on  a  mass  basis)  of  the 
gas  mixture  presented  in  reference  15  at  temperatures  near  the  middle  of  the  range  of 
validity  presented  in  figure  3.  At  the  low  and  high  temperature  extremes  of  the  range  of 
validity,  the  respective  percentages  of  this  gas  mixture  accounted  for  are  approximately 
80  and  70.  Hence,  the  presently  assumed  pyrolysis  gases  are  considered  to  be  a  reason¬ 
able  approximation  for  the  case  where  the  cracking  reactions  approach  equilibrium. 

If,  on  the  other  hand,  the  cracking  reactions  were  kinetically  slow  so  that  large  molecular 
weight  hydrocarbon  species  reached  the  char  surface,  the  pyrolysis  gases  in  assumption 
(f)  would  be  incorrect  and  the  presence  of  heavy  gases  would  cause  the  boundary -layer 
gas  mixture  to  deviate  still  further  from  an  effectively  binary  mixture.  The  fact  that  the 
major  pyrolysis  gas  species  measured  experimentally  (refs.  13  and  14)  and  calculated 
(ref.  15)  are  the  same  strongly  suggests  that  char  layer  cracking  reactions  approach 
equilibrium. 


Physical  Interpretation  of  Approximate  Solution 

A  physical  interpretation  of  the  approximate  solution  may  be  obtained  from  an 
examination  of  equations  (128)  and  (129).  The  numerator  of  equation  (128)  is  the  well- 
known  expression  for  the  diffusion-controlled  oxygen  flux  and  is  the  maximum  rate  at 
which  oxygen  can  be  transported  across  the  boundary  layer  to  the  reacting  char  surface. 
Since  equation  (128)  is  an  expression  for  the  steady-state  ablation  rate,  fW*  must  rep¬ 
resent  the  mass  of  oxygen  used  in  ablating  a  unit  mass  of  virgin  plastic.  If  every  carbon 
atom  initially  present  in  the  virgin  plastic  is  assumed  to  react  with  the  oxygen  in  the  vir¬ 
gin  plastic  to  form  CO,  then  the  term  - v  -  cp -  represents  the  number 

crM^  +  tMh  +  i^Mn  +  <pMQ 

of  moles  of  carbon  per  unit  mass  of  virgin  plastic  available  to  react  with  oxygen  diffusing 
through  the  boundary  layer.  Hence  fW*  (eq.  (129))  represents  the  mass  of  oxygen  from 
external  sources  required  to  react  all  the  carbon  in  a  unit  mass  of  virgin  plastic  to  pro¬ 
duce  CO.  This,  then,  is  the  physical  interpretation  of  the  approximate  solution:  all  the 
carbon  present  in  the  virgin  plastic  is  reacted  to  form  CO  and  the  reaction  rate  is  deter¬ 
mined  by  the  rate  at  which  oxygen  can  be  supplied  by  pyrolysis  of  the  virgin  plastic  plus 
diffusion  through  the  boundary  layer.  It  is  in  this  respect  that  the  present  analysis  differs 
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most  significantly  from  previous  analyses  of  char  oxidation  and  carbon  combustion.  In 
these  previous  analyses,  only  the  carbon  present  in  the  char  layer  is  considered  reactable. 

Effect  of  Virgin  Plastic  Composition  on  Boundary -Layer  Composition 
at  the  Char  Surface,  Char  Oxidation  Rate,  Ablation  Rate, 
and  Effective  Heat  of  Ablation 

The  effect  of  virgin  plastic  composition  on  pyrolysis  rate  is  shown  by  equa¬ 
tions  (128)  and  (129).  However,  the  transfer  coefficient  h  which  appears  in  equa¬ 
tion  (128)  is  dependent  on  the  virgin  plastic  hydrogen  content  as  shown  in  equation  (135). 
When  equation  (135)  is  substituted  into  equation  (128),  the  following  relation  is  obtained: 


fW%Ko>e„ 


£W*+  0.7(1  +lSmH/  Kq  e 


(139) 


Equation  (139)  indicates  that  mvp  decreases  with  increasing  fW*  and  increasing 
transpiration  factor  rj.  Since  fW*  increases  with  increasing  virgin  plastic  carbon 
content  and  decreasing  oxygen,  nitrogen,  and  hydrogen  content  and  r]  increases  with 
increasing  virgin  plastic  hydrogen  content,  the  question  arises  as  to  whether  the  lowest 
pyrolysis  rate  would  be  obtained  from  a  material  with  a  high  carbon  content  or  a  high 
hydrogen  content.  In  an  attempt  to  answer  this  question,  equation  (139)  was  used  to  make 
calculations  of  mVp  for  a  material  which  contained  only  carbon  and  hydrogen.  (Eq.  (129) 
shows  that  the  presence 


of  oxygen  or  nitrogen  in 
the  virgin  plastic  always 
produces  higher  pyroly¬ 
sis  rates.)  The  results 
of  these  calculations  are 
presented  in  figure  4  as 
plots  of  mVp/h° 
against  m^  for 
kQ  e  =  0.050  and  0.232 
and  l  =  1/3.  For  a 
specified  flow  environ¬ 
ment  ^h°  and  Kq  e^, 
the  minimum  value  of 
mVp  is  obtained  from 
a  material  which  is 


Carbon  mass  fraction  of  virgin  plastic  .m^ 

Figure  4.-  Variation  of  dimensionless  steady-state  ablation  rate  with  carbon  mass  fraction  for 
plastic  which  contains  only  carbon  and  hydrogen. 
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composed  entirely  of  carbon.  Numerous  investigators  have  impregnated  porous  carbon 
materials  with  hydrocarbon  resins  in  the  hope  that  the  transpiration  cooling  effect  pro¬ 
duced  by  the  vaporization  of  the  resin  would  reduce  the  mass  loss  rate  of  the  material. 
The  results  of  this  study  show  that,  within  the  framework  of  the  assumption  presented 
earlier,  the  effect  of  impregnating  a  porous  carbon  with  hydrocarbon  resins  would  be  to 
increase  rather  than  reduce  the  mass  loss  rate. 

In  order  to  illustrate  the  effects  of  virgin  plastic  composition,  calculations  of  char 
oxidation  rate,  ablation  rate,  and  char  surface  boundary -layer  composition  were  carried 
out  for  four  typical  hydrocarbon  resins  and  graphite.  The  empirical  formulas  and 
associated  parameters  used  in  these  calculations  are  presented  in  table  I.  The  indices 
in  the  empirical  formulas  have  been  adjusted  to  give  the  numbers  of  carbon  and  hydrogen 
atoms  per  atom  of  oxygen  and  the  materials  are  listed  in  order  of  decreasing  carbon  con¬ 
tent.  As  pointed  out  previously  the  parameter  fW*  represents  the  mass  of  oxygen 
required  to  ablate  a  unit  mass  of  virgin  plastic.  (For  graphite,  no  pyrolysis  is  involved 
and  hence  fW*  represents  the  mass  of  oxygen  required  to  oxidize  a  unit  mass  of 
graphite.)  The  values  presented  in  table  I  show  that  fW*  increases  markedly  as  the 
relative  carbon  content  of  the  material  increases.  The  symbol  mH  represents  the 
mass  of  hydrogen  per  unit  mass  of  virgin  plastic  and  is  related  to  the  transpiration 
cooling  factor  t]  by  equation  (125).  The  significant  point  to  be  made  with  the  values  pre¬ 
sented  in  table  I  is  that  while  mH  differs  significantly  for  the  five  materials  considered, 

TABLE  I 

EMPIRICAL  FORMULAS  AND  ASSOCIATED  PARAMETERS  FOR  GRAPHITE 
AND  FOUR  TYPICAL  HYDROCARBON  PLASTICS 


Material 

Empirical 

formula 

fW* 

mH 

V 

(a) 

Graphite 

C 

1.333 

0 

0.700 

Epoxy 

resin 

C6.5  H7.34  °1 

.869 

.0725 

.874 

Phenolic 

resin 

C5.16  H5  °1 

.802 

.0602 

.849 

Polyester 

resin 

C2.5  H2  °1 

.500 

.0416 

.810 

Cellulose 

C1  2  H,  Oj 

.099 

.0617 

.851 

aIn  the  calculation  of  7]  (see  eq.  (125)),  the  exponent  l  was  assumed 
to  be  1/3. 
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the  maximum  variation  in  77  is  only  25  percent.  This  again  indicates  that  large 
increases  in  the  transpiration  cooling  effect  will  not  be  obtained  by  increasing  the  hydro¬ 
gen  content  of  a  plastic. 

By  using  the  empirical  formulas  presented  in  table  I  and  equations  (130),  (131), 
and  (132),  the  boundary -layer  compositions  at  the  char  surface  during  the  steady  abla¬ 
tion  of  cellulose,  polyester,  phenolic,  and  epoxy  resins,  and  a  completely  carbonaceous 
material  were  calculated  and  are  presented  in  figure  5  as  variations  of  N9,  H9,  or  CO 


mass  fraction  with  free -stream  oxygen 
mass  fraction.  The  char  surface 
boundary -layer  compositions  are  com¬ 
parable  for  all  the  materials  except 
cellulose.  Cellulose,  with  its  high  oxygen 
and  hydrogen  contents,  produces  much 
larger  quantities  of  CO  and  than  the 
other  materials  do.  The  curves  pre¬ 
sented  in  figure  5  show  that  (with  the 
exception  of  cellulose)  the  boundary - 
layer  composition  for  the  plastics  is 
similar  to  that  produced  by  graphite 

i-e-’  KH2  <<  Kco  and  kn2)-  °ne  °f 

the  principal  assumptions  involved  in  the 
present  analysis  is  that  the  unity  Lewis 
number  heat-mass  transfer  analogy  is 
valid.  It  is  generally  agreed  that  the 
unity  Lewis  number  assumption  is  rea¬ 
sonably  valid  for  carbon  oxidation. 

(See  ref.  9.)  The  comparable  boundary - 
layer  composition  shown  in  figure  5  for 
epoxy  and  phenolic  resins  and  graphite 
suggests  that  the  unity  Lewis  number 
assumption  may  be  a  reasonable  approx¬ 
imation  for  the  present  case. 


Graphite  C(s) 


Free-stream  oxygen  mass  fraction,  Kq  0 

Figure  5.-  Effect  of  virgin  plastic  composition  on  the  variation  of 
boundary- layer  composition  at  the  char  surface  with  free-stream 
oxygen  mass  fraction. 


The  effect  of  virgin  plastic  composition  on  steady-state  ablation  rate  is  illustrated 
in  figure  6,  where  ratios  of  ablation  rate  to  Kn  h°  are  plotted  against  oxygen  mass 
fraction,  Kq  0.  The  data  are  presented  in  this  manner  since  the  product  Kq  0h  rep¬ 
resents  the  maximum  rate  at  which  oxygen  can  be  transported  by  diffusion  to  the  char 

*  o 

surface  and  hence  characterizes  the  flow  environment.  The  ratio  gh  indi¬ 

cates  the  ablator's  vulnerability  to  oxidation.  The  curve  for  graphite  is  included  since  it 
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represents  (as  demonstrated  earlier)  the  min¬ 
imum  ablation  rate  attainable  by  a  hydrocarbon 
material.  The  curves  presented  in  figure  6 
show  that  plastics  exhibit  a  range  of  carbon 
content  which  causes  significant  differences  in 
char  oxidation  rates  and  ablation  rates.  The 
ablation  rates  calculated  for  phenolic  and 
epoxy  resins  (the  two  plastics  most  commonly 
used  in  ablative  heat  shields)  are  about  1.4  to 
1.6  times  as  large  as  the  theoretical  minimum 
rates  (for  graphite),  and  the  values  for  poly¬ 
ester  and  cellulose  exceed  those  for  graphite 
by  factors  of  from  2.7  to  3  and  5  to  12, 
respectively.  From  time  to  time  it  has  been 
suggested,  in  the  literature,  that  wood  (which 
is  primarily  cellulose)  might  be  an  acceptable 
ablator.  The  calculations  presented  in  fig¬ 
ure  6,  however,  indicate  that  cellulose  would 
not  be  an  acceptable  material  for  flight 
regimes  where  steady  ablation  and  diffusion- 
controlled  char  oxidation  are  obtained.  Of 
the  three  resins,  phenolic,  epoxy,  and  polyester,  the  present  results  indicate  that  phenolic 
and  epoxy  would  be  the  least  vulnerable  to  oxidation  and  indeed  these  are  the  two  resins 
most  commonly  used  in  ablative  heat  shields.  Even  though  phenolic s  and  epoxys  are 
relatively  efficient  ablators,  the  curves  presented  in  figure  6  predict  that  completely 
carbonaceous  materials  (such  as  prepyrolyzed  resins)  would  have  mass-loss  rates  that 
were  from  60  to  70  percent  of  the  rates  obtainable  from  phenolics  and  epoxys.  Most 
carbonaceous  materials  have  relatively  high  thermal  conductivities  and  densities  and 
hence  are  poor  insulators.  The  proper  interpretation  of  the  curves  presented  in  fig¬ 
ure  6  is,  then,  that  if  completely  carbonaceous  materials  with  insulating  capabilities 
comparable  to  those  of  plastics  could  be  developed,  then  these  carbonaceous  materials 
would  be  superior  ablators  by  virtue  of  these  lower  mass  loss  rates. 

Estimates  of  the  weight  of  material  ablated  by  a  given  aerodynamic  heat  load  can  be 
obtained  by  defining  an  "effective  heat  of  ablation."  For  present  purposes  a  cold-wall 
effective  heat  of  ablation  is  used  and  is  defined  as  follows: 
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Figure  6.-  Effect  of  virgin  plastic  composition  on  the  ratio  of 
ablation  rate  (or  char  oxidation  rate)  to  diffusion  oxygen 
flux. 
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Substitution  of  equation  (139)  into  equation  (140)  yields 


Qeff  ” 


iH. 


(141) 


For  diffusion-controlled  oxidation,  the  response  of  the  ablator  is  determined  by  the  flux 
of  oxygen  to  the  char  surface  and  is  completely  independent  of  the  aerodynamic  heat  flux. 
As  a  result,  equation  (141)  contains  no  thermal  properties  such  as  specific  heats,  thermal 
conductivities,  or  heats  of  pyrolysis.  It  is  only  in  mass  transport  controlled  cases,  such 
as  the  present  one,  that  an  effective  heat  of  ablation  such  as  that  defined  in  equation  (140) 
has  meaning.  Where  the  response  of  the  ablator  is  determined  by  aerodynamic  heat  flux, 
a  thermochemical  heat  of  ablation  such  as 


m 


vp 


(142) 


^  where  qne^  is  the  heat  flux  to  a  nonablating  calorimeter  that  has  the  same  tempera 
ture,  emissivity,  and  catalytic  efficiency  as  the  ablator^  must  be  used. 


Equation  (141)  was  used  to  calculate  values  of  Qeff  for  each  of  the  materials 
appearing  in  table  I  for  KQ  =  0.232,  which  corresponds  to  air.  The  variation  of  these 
effective  heats  of  ablation  with  stagnation  enthalpy  is  presented  in  figure  7,  which  shows 


_ — —  Completely  carbonaceous  (graphite) 

Btu  - Ep°*y 


lb  - Phenolic 


Figure  7.-  Effect  of  virgin  plastic  composition  on  cold-wall  effective  heat  of  ablation. 


35 


that  the  cold-wall  effective  heats  of  ablation  of  epoxy,  phenolic,  polyester,  and  cellulose 
are,  respectively,  72,  67,  45,  and  19  percent  of  that  for  a  completely  carbonaceous  mate¬ 
rial.  Since,  for  a  given  cold-wall  heat  load,  the  weight  of  material  ablated  is  inversely 
proportional  to  He^,  figure  7  may  be  interpreted  as  showing  that  the  respective  weights 
of  material  ablated  from  heat  shields  composed  of  epoxy,  phenolic,  polyester,  and  cellu¬ 
lose  would  be  1.39,  1.49,  2.22,  and  5.25  times  the  weight  of  material  ablated  from  a  com¬ 
pletely  carbonaceous  shield. 

Effect  of  Pyrolysis  Gas  Composition  on  the 
Heat  Release  Due  to  Chemical  Reaction 

In  order  to  illustrate  the  effects  of  varying  pyrolysis  gas  composition,  calculations 
were  made  for  two  assumed  pyrolysis  gas  compositions.  These  compositions,  which  are 
presented  in  table  II,  were  based  on  data  presented  in  references  13  and  14.  The  results 

TABLE  II 

ASSUMED  PYROLYSIS  GAS  COMPOSITIONS 


Item 

Composition  1 
(phenolic  resin) 

Composition  2 
(phenolic  nylon) 

XR  . 

H2 

0.5199 

0.2061 

o 

o 

IX 

0.1954 

0.2875 

. 

2 

0.0013 

0 

*co2  . 

0 

0.0222 

XCH4  . 

0.1427 

0.1980 

xc2h2 . 

0 

0.1340 

XC2H4 . 

0 

0.0006 

xnh3  . 

0.1407 

0.1515 

xh2o . 

0 

0 

fW* . 

0.724 

0.601 

mH . 

0.0974 

0.0795 

1 . 

0.921 

0.886 

CqHjN^O^  .  .  . 

c5.9  h10.37  n0.73  °1 

c3.74  h5.80  n0.46  °1 

f . 

0.462 

0.23 
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of  these  calculations  are  compared  with  those  for  graphite  and  are  presented  in  fig¬ 
ures  8,  9,  and  10. 

From  table  II  it  is  seen  that  the  two  compositions  differ  significantly  with  respect 
to  H2,  CO,  and  C2H2  content.  A  comparison  of  table  H  with  table  I  shows  that  the  values 
of  fW*  and  the  virgin  plastic  carbon  and  oxygen  contents  derived  from  the  pyrolysis  gas 
compositions  are  lower  than  those  derived  from  the  molecular  structure  of  phenolic 
resin.  Because  of  the  extreme  difficulty  of  experimentally  determining  pyrolysis  gas 
composition,  and  because,  in  the  compositions  presented  in  table  II  trace  species  and 
condensibles  included  in  references  13  and  14  were  neglected,  the  compositions  presented 
in  table  II  must  be  considered  as  assumed  rather  than  experimentally  measured.  Con¬ 


sidered  in  this  light,  the  agreement  in  fW 

Figure  8  presents  the  boundary - 
layer  composition  at  char  surface  as  a 
function  of  Kq  0  for  two  assumed  com¬ 
positions  and  graphite.  Comparing  fig¬ 
ure  8  with  figure  5  shows  that  the  char 
surface  boundary-layer  composition  pro¬ 
duced  by  the  assumed  pyrolysis  gases 
closely  approximates  that  derived  from 
the  molecular  structure  of  phenolic 
resin. 

In  figure  9  the  variation  of  the 
maximum  heat  of  reaction  (per  pound  of 
char)  with  char  surface -temperature  is 
presented  for  graphite  and  the  two  pyrol¬ 
ysis  gas  compositions.  The  maximum 
heat  of  reaction  AHR  for  composition  1 
is  only  slightly  lower  than  that  for  graph¬ 
ite,  whereas  AHR  for  composition  2  is 
approximately  twice  that  for  graphite.  A 
large  part  of  this  difference  is  due  to  the 
difference  in  f  (mass  of  char  produced 
when  a  unit  mass  of  virgin  material  is 
pyrolyzed)  between  compositions  1  and  2. 
(See  eq.  (116).)  Figure  9  shows  that  AH^ 
(per  pound  of  char)  can  be  a  strong  function 
of  virgin  plastic  composition  and  pyrolysis 
gas  composition. 


is  perhaps  not  too  bad. 


o 

0) 

CL 


Figure  8.-  Variation  of  ooundary- layer  composition  at  char  surface 
with  free-stream  oxygen  mass  fraction  as  calculated  from 
assumed  pyrolysis  gas  compositions. 
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□  Phenolic  resin  (composition  I  ) 

O  Phenolic  nylon  (composition  2  } 


Char  surface  temperature  ,  Tw 

Figure  9.-  Variation  of  maximum  heat  of  reaction  with  char  surface  temperature  as  calculated  from  assumed 

pyrolysis  gas  compositions. 


As  shown  by  equations  (117),  (118),  (119),  and  (134)  an  effective  enthalpy  potential 
due  to  chemical  reaction  can  be  defined  as  a  function  of  virgin  plastic  composition,  pyrol¬ 
ysis  gas  composition,  and  free-stream  oxygen  mass  fraction.  This  effective  enthalpy 
potential  ^per  unit  KQ  ^  is  plotted  against  char  surface  temperature  in  figure  10.  Note 


J 

kg  Bfu 


IOOO  2000  3000  4-000  5000  6000  7000  °R 


Char  surface  temperature  ,  Tw 

Figure  10.-  Variation  of  effective  enthalpy  potential  due  to  oxidation  with  char  surface  temperature  as 
calculated  from  assumed  pyrolysis  gas  compositions. 
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that  AHch  for  composition  2  is  larger  than  that  for  graphite  whereas  AHch  for  com¬ 
position  1  is  less  than  that  for  graphite.  The  parameters  responsible  for  these  differ¬ 
ences  can  be  identified  by  examining  equation  (134).  Rearrangement  of  equation  (134) 
yields 


iHCh  = 


(4Hch) 


ej 

C(s)  M()  W  -  «: 


X 


f,C0  sco  -  25co2  +  sh2o 


(143) 


where  the  first  term  on  the  right-hand  side  is  the  effective  enthalpy  potential  for  the  com 
bustion  of  graphite.  Note  that  the  second  term  on  the  right-hand  side  will  be  positive  or 
negative  ^and  hence  AH^  will  be  greater  than  or  less  than  ^AH^^  j  depending  on 


X  ^,1 


whether  AH^  is  greater  than  or  less  than  = 


Also  note  that 


xco  +  2xco2  +  xh2° 


for  a  given  value  of  AH^  qq  - 


xco  +  2xco2  +  xh2o 


the  magnitude  of  AH  ^  depends 


on  the  size  of  the  ratio  ^  Since  AH^  qq  is  a  negative  quantity,  AH^  will  be 


larger  than  ^AHC^  ^  when  ^  X^AH^  i  is  Positive.  Of  all  the  pyrolysis  gas  con¬ 
stituents  considered  in  this  paper,  only  C2H2  and  CgH^  have  positive  values  of  AH. 

Hence  AH^  becomes  larger  as  the  amounts  of  and  CgH^  in  the  pyrolysis  gases 

increase.  Note  from  figure  10  and  table  II  that  AH^  for  pyrolysis  gas  composition  2 
^  which  contains  larSer  than  AH^  ^or  Pyro^ys^s  Sas  composition  1  ^  which 

does  not  contain  C2H2).  Vojvodich  and  Pope  reported  in  reference  17  heating  due  to 
chemical  reactions  much  larger  than  could  be  accounted  for  by  carbon  combustion.  From 
equation  (143)  it  is  seen  that  very  large  values  of  AH^  could  be  produced  by  a  plastic 
which  had  a  small  value  of  (a  -  cp)  and  produced  pyrolysis  gases  containing  large  amounts 
of  C2H2  and  C2H^. 


39 


SUMMARY  OF  RESULTS 


An  analytical  study  of  diffusion-controlled  char  oxidation  and  its  effect  on  the  abla¬ 
tion  of  hydrocarbon  plastics  has  been  carried  out.  Although  the  solution  obtained  is 
based  on  a  number  of  simplifying  assumptions,  it  is  believed  to  present  an  accurate  con¬ 
ceptual  picture  of  the  diffusion -controlled  char -oxidation  process.  The  steady-state 
expressions  resulting  from  this  analysis  can  be  used  directly  to  provide  a  clear  physical 
picture  of  char  oxidation.  This  analysis  has  also  produced  general  expressions  which 
can  be  used  in  conjunction  with  a  transient  heat  conduction  computer  program  to  yield 
transient  ablation  and  char  oxidation  rates,  but  such  a  procedure  does  not  lend  itself  to 
demonstrating  the  general  characteristics  of  diffusion -controlled  char  oxidation.  Accord¬ 
ingly,  the  discussion  and  results  presented  in  the  present  report  were  restricted  to  the 
steady-state  case.  The  results  of  this  study  are  summarized  as  follows: 

1.  The  present  analysis  differed  from  previous  analyses  of  diffusion -controlled 
char  oxidation  in  one  important  respect:  previous  investigators  have  assumed  that  only 
the  carbon  in  the  char  layer  would  react  and  that  the  rate  of  reaction  would  be  determined 
by  the  diffusion  rate  of  oxygen  through  the  boundary  layer  but  the  present  analysis  has 
shown  that  the  carbon  present  in  both  the  char  and  the  pyrolysis  gases  (i.e.,  all  the 
carbon  originally  present  in  the  virgin  plastic)  reacts  and  that  the  reaction  rate  is 
determined  by  the  rate  at  which  oxygen  can  be  supplied  by  diffusion  through  the  boundary 
layer  plus  pyrolysis  of  the  virgin  plastic.  Because  of  this  fact,  the  char  oxidation  rate, 
the  steady-state  ablation  rate,  the  boundary-layer  composition  at  the  char  surface,  and 
the  heat  released  by  chemical  reactions  were  found  to  be  dependent  on  the  composition 

of  the  virgin  plastic  and  the  pyrolysis  gases. 

2.  For  a  given  flow  environment  (i.e.,  oxygen  mass  fraction  and  heat-transfer  coef¬ 
ficient)  the  char  oxidation  rate  and  the  pyrolysis  rate  were  found  to  increase  as  the 
hydrogen,  nitrogen,  and  oxygen  contents  of  the  virgin  plastic  increased  and  to  decrease 
as  the  carbon  content  of  the  virgin  plastic  increased.  Even  when  transpiration  cooling 
effects  were  accounted  for,  the  lowest  steady-state  pyrolysis  rate  was  obtained  from  a 
material  composed  entirely  of  carbon. 

3.  The  cold-wall  effective  heat  of  ablation  of  a  plastic  material  was  found  to  be 
determined  by  the  free -stream  oxygen  mass  fraction,  the  stream  enthalpy  potential,  and 
relative  number  of  carbon,  hydrogen,  nitrogen,  and  oxygen  atoms  in  the  virgin  plastic. 
Example  calculations  carried  out  for  a  completely  carbonaceous  material,  epoxy,  phenolic 
and  polyester  resins,  and  cellulose  showed  that  the  effective  heats  of  ablation  for  the 
polymers  would  be,  respectively,  72,  67,  45,  and  19  percent  of  the  effective  heat  of  abla¬ 
tion  calculated  for  the  carbonaceous  material. 
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4.  The  maximum  heat  of  reaction  per  pound  of  char  was  a  strong  function  of 
pyrolysis  gas  composition  and  the  fraction  of  virgin  material  which  was  converted  to 
gas,  and  the  heat  released  by  both  gas  phase  and  heterogeneous  reactions  can  be 
accounted  for  by  an  "effective  enthalpy  potential"  which  is  a  function  of  the  stream 
oxygen  mass  fraction,  the  relative  number  of  carbon,  oxygen,  hydrogen,  and  nitrogen 
atoms  in  the  virgin  plastic  and  the  composition  of  the  pyrolysis  gases.  The  pyrolysis 
gases  considered  in  this  investigation  were  CO,  C02,  H2,  N2,  NH^,  H20,  CH^,  CgHg, 
and  C2H4.  Both  the  maximum  heat  of  reaction  and  the  effective  enthalpy  potential  were 
found  to  increase  with  increasing  C2H2  and  C2H4  content.  Particularly  high  heat 
releases  were  predicted  for  polymers  which  contained  nearly  equal  numbers  of  carbon 
and  oxygen  atoms  and  produced  pyrolysis  gases  containing  large  quantities  of  C2H2  and 
and  C2H4- 
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